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ABSTRACT 

We present the first scientific results from ttie Sydney- A AO Multi-Object IFS (SAMI) at the Anglo- 
Australian Telescope. This unique instrument deploys 13 fused fibre bundles (hexabundles) across a 
one-degree field of view allowing simultaneous spatially-resolved spectroscopy of 13 galaxies. During 
the first SAMI commissioning run, targeting a single galaxy field, one object (ESO 185-G031) was 
found to have extended minor axis emission with ionisation and kinematic properties consistent with a 
large-scale galactic wind. The importance of this result is two-fold: (i) fibre bundle spectrographs are 
able to identify low-surface brightness emission arising from extranuclear activity; (ii) such activity 
may be more common than presently assumed because conventional multi-object spectrographs use 
single-aperture fibres and spectra from these are nearly always dominated by nuclear emission. These 
early results demonstrate the extraordinary potential of multi-object hexabundle spectroscopy in 
future galaxy surveys. 
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1. INTRODUCTION 

1.1. Feedback and Galactic Winds 

Feedback, from stars and active galactic nuclei (AGN), 
plays a vital role in regulating star formation in galaxies, 
thus greatly impacting their growth and evolution. One 
of the most spectacular feedback phenomena observed 
throughout the cosmos are galactic winds. Their effects 
can be seen across the electromagnetic spect rum from 



radio waves (jDuricet al. 1983), mid-infrared (Roussel et 
'al.||2QlQ[), and UV ([Hoopes et al. 2005) through to X 



et al. 



raysH ^nckland et al. ||2004| , and even gamma rays (Su 



2010| ). The importance of winds in galactic evom 



tion IS still unclear although they are suspecte d to play 
a rol e in the observed mass-met allicity relation (Dekel & 
Silk ll986) and the reduced baryon fra ction in galaxies 



comp ared to the universal background ( [McGaugh et al. 
2010). Outflows are routinely invoked in numerical sim- 



ulations to forcibly remove the build-up of low angular 
momentum material at the centre of syn thetic galaxies 
dBe rtone et"ar||2005[ [Brook et al. | |2010[ [Broo ks 2010). 
At the present time, this appears to be the only viable 



mechanis m for producing exponential disks ([Sharma et 
ar|[2QT2l ) 
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Large-scale winds are common in starburst and active 
galaxies and are therefore likely to be important at high 
redshifts where starbursts dominate sta r formation and 
AGN activity was at its peak ( Veilleux et al. ||2QQ5 ). In 
general, winds are highly energetic and largely mvisible, 
where their presence is betrayed by material entrained 
from the galactic disk. This is evident from the slow ro- 
tation of ou tflowing fllaments about the wind axis ( [Shop-] 
beh & Bland-Hawthorn 1998; Greve 2004). 



I'here has never been a systematic survey to assess 
the statistical occurrence of galactic wind sources in any 
waveband. While integrated spectra can show evidence 
of outflows, very detailed observations are required to 
properly separate the f aint wind material from the back- 
ground galaxy ( [Cecil et al. 2001). The clearest indica- 
tion of a wind typically comes from either (i) s patially 
resolved X-ray imaging (Strickland et al. 2004); or (ii) 
combining spatially re solved ionization diagnost ics and 
kinematic in formation (jVeilleux & Rupke 2002; Sharp & 
Bland-Ha wthorn , as we demonstrate here. 

Because of the need for detailed spatially-resolved 
data, it has been diflicult to envisage how a large-scale 
systematic survey of galactic winds would be possible 



if t hese a mount to only 1% of all galaxies ( [Veilleux et 



ah] 2005[). With an instrument concept such as the 
wd: 



Sydney- A AO Multi-ob ject Integral field spectrograph 
(IFS) (SAMI; [Groom e t al. 2012) a survey of this kind 
is now possible. As argued below, we envisage surveys of 
thousands of galaxies, where every object has spatially 
resolved IFS data, enabling the detection of extended, 
low surface brightness material and thus the identifica- 
tion of complex structures and spatial gradients in galaxy 
properties. 

1.2. The Sydney-AAO Multi-object Integral field 
spectrograph (SAMI) 
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The Sydney- A AO Multi-object IFS is an innova- 
tive new instrument on the Anglo- AustraUan Telescope 
(AAT), providing 13 deployable fused fibre bundles, 
called hexabundles, over a large field of view. One of 
the strongest science motivations behind such an instru- 
ment is the study of complexity in galaxies. Understand- 
ing complex processes, in particular the fiow of gas in 
and out of galaxies, is vital for solving the major prob- 
lems of galaxy formation and evolution. Phenomena 
such as accretion, feedback, and the dynamical distur- 
bance of gas through external and secular processes can 
manifest themselves in subtle ways via spatial gradients 
and/or structure in chemical abundances, ionisation, and 
kinematic properties, to name a few. Finding such sub- 
structure in galaxy properties is not trivial and relies on 
spatially-resolved spectral data. It is generally not pos- 
sible for single-aperture spectroscopic studies to identify 
what kinds of objects exhibit these features. The dis- 
covery space probed with multiplexed IFS allows us to 
simultaneously find and study such objects in detail. 

SAMI is the first large-scale multiplexed IFS to use 
hexabundles, a new optical fibre imaging bundle de- 
signed an d devel oped at the University of Sydney (Bland- 
[Hawthorn et al. ||2Qll| [Bryant et al. |2Qll[ ). Hexabundles 
consist ot optical tibres with reduced cladding arranged 
in a circular pattern and lightly fused over a short length, 
yielding an imaging bundle wi th a high fill factor and ex- 
cellent optical performance (Bryant et al. |2011). SAMI 
uses 13 hexabundles deployable across a 1-degree field of 
view with each of the SAMI bundles containing 61 fibres. 
The instrument also contains 26 separate single fibres for 
sky subtractio n. SAMI feeds the double-beam AAOmega 
spectrograph ( [Sharp et al. [2006 ). 

SAMI is thus designed to simultaneously observe 13 
spatially extended objects (e.g. galaxies), with 61 spa- 
tial samples across each target and spectral information 
for each of these spatial samples. The power of this tech- 
nique is immediately apparent if we consider two things - 
the aperture effects inherent in single fibre galaxy surveys 
and the time-consuming nature of current IFS surveys. 

The effects of a single fixed aperture size (e.g. a fibre on 
sky) on large galaxy surveys have been known for some 
time. The problem is simple: the properties derived from 
aperture spectra of galaxies depend on both the position 
and size of the aperture one is observing with. For galax- 
ies larger than the aperture only the central part of the 
galaxy is observed and some of the light in the outskirts 
of the object is lost. Since the apparent size of a galaxy 
on sky will vary with redshift, even over modest redshift 
ranges, this makes it very difficult to compare samples 
across different epochs and could affect the interpretation 
of the evolution of global galaxy properties. Likewise, if 
the observing aperture is not well centred on a galaxy it 
is difficult to say what relevance the derived properties 
will have to the object's true integrated properties. This 
is particularly relevant if a galaxy is not morphologically 
well-behaved, for example, how is one to determine the 
"centre" position of a merging or otherwise disturbed 
system? 

This problem is addressed i n various ways using com- 
plex aperture corrections (e.g. Brinchmann et al. (2004) 
for SDSS). While these aperture corrections can be 
shown to be broadly correct when applied to a popula- 
tion of galaxies (albeit with significant scatter remaining 



Date 


Run # 


Total Exposure Time (s) 


Format 


2nd July 


29-31 


7200 


Pointed 


3rd July 


29-37 


14034 


Dithered 



Table 1 

SAMI observations of the wind galaxy. 



in the derived properties, see Hopkins et al. (2003)), they 
can suffer gross inaccuraci es when dealing ^ with i ndivid- 



ual galaxies. For example Gerssen et al. (2012) found 



that for a sample of star- form ing galaxies in SPSS, the 



star formation rates derived by Brinchmann et al. (2004) 
were underestimated by a median factor of 2.5, with a 
scatter of 175%. As well as this need for aperture cor- 
rection, single fibre spectra are inherently limited as they 
can only reveal the average properties of galaxies. In con- 
trast to IFS data, any spatial variation across the object 
itself is lost. 

Although the need for them is clear, large surveys of 
galaxies with IFS data have heretofore been difficult and 
time-consuming to execute. This is because most IFS are 
single-object instruments, meaning only one galaxy can 
be observed at a time. The largest sample to date is that 
of the ATLAS-3D project ( Capp ellari et al. |201lD , num- 
bering 260 early-type galaxies each observed with the 



et al 



SAURON spectrograph. The CALIFA project ( Sanchez 
2012!) aims to build a sample of 600 galaxies usmg 



the PMAS IFS. A t higher redshift the IMAGES survey 
( Yang et al. |[2QQ8| used the FLAMES/GIRAFFE MOS- 
Il^'U system to observe a sample of 63 galaxies at reds hifts 
z^O. 4-0^ 75, while the SINS survey ( Forster-Schreiber et| 
ar] |2QQ9| used the SINFONI IFS to observe a sample 
of d3 galaxies at redshifts z~1.2-3.6. These samples 
clearly demonstrate the richness and potential of large 
IFS samples but were extremely time-consuming to ob- 
serve. With its multiplex factor of 13, SAMI will be able 
to observe a sample of ~350 galaxies in twenty nights. 
The first SAMI com missioning results were published 



Groom et al. (2012). To further demonstrate the sci 
entitle capability of SAMI, here we take a closer look at 
just one of the objects observed during commissioning 
in July 2011. Specifically we present observations of a 
serendipitously discovered galactic wind, whose detection 
was enabled by spatially resolved SAMI spectroscopy of 
the galaxy in question. This discovery demonstrates the 
power of multi-object IFS to identify rare objects that 
are impossible to detect in single integrated spectra. 

Section ^describes our instrument set-up and observa- 
tions, while Section [3] describes our data reduction proce- 
dure with particular attention paid to the dithered SAMI 
observation. Section [4] describes our results, with Section 
[5] laying out our conclusions. Ou r analysis uses the cos- 
mological parameters provided by Komatsu et al. (2009) 
with (^A=0.723, (^M=0.277 and Ho=70.2 kms-'Mpc^ 

2. OBSERVATIONS 

SAMI was mounted at the prime focus of the AAT 
behind the triplet corrector which provides a 1-degree 
field of view with a plate scale of 15.2'Ymm. Each of the 
SAMI hexabundles contains 61 105/im-diameter fibres in 
a circular pattern, meaning each hexabundle has a field 
of view of 14^9 on sky, with a sampling of V!6 per fibre. 
From its position at prime focus, a 42m fibre cable runs 
from SAMI to the AAOmega spectrograph in the coude 
room. 
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ESO 185-G031 



Figure 1. An R-band image of ESO 185-G031 taken with the 
Maryland-Magellan Tunable Filter (MMTF) instrument on the 6.5- 
meter Walter Baade Magellan telescope at Las Campanas Obser- 
vatory. The galaxy is a gas-rich spiral with well-defined arms. The 
near-side arm is clearly visible running south-east to north-west 
along the disk. The red circle indicates the footprint of the SAMI 
hexabundle on the galaxy. 

AAOmega is a double-beam spectrograph which can 
operate with a range of dichroics and gratings, yielding 
a corresponding range of values for resolution and wave- 
length coverage in each spectrograph arm. For these ob- 
servations we used the 580V grating in the blue arm, 
providing spectral coverage from 3700-5700 A and a res- 
olution of = A/ A A 1730. This grating was chosen to 
cover the most interesting spectral features in the blue 
end of the spectrum, including the [O ll] A 3727 A emis- 
sion feature and the "4000- Angstrom" continuum break, 
D4000- However, since the type of fibre used for the 42m- 
long SAMI fibre run was shown to have poor throughput 
in the blue end of th e spectrum (this fibr e run is currently 
being upgraded, see Groom et al. (2012j) for more details) 
we did not observe any usable data below 4000 A, with 
the throughput beginning to decline for A < 4500 A. We 
did, however, robustly measure other features redward 
of 4500 A, including, for this galaxy, both the H/3 and 
[Oiii] AA4959, 5007 A doublet emission features. 

In the red arm, we used the lOOOR grating provid- 
ing spectral coverage over the range 6250-7350 A at a 
resolution of ~ 4500. This grating was chosen to 
provide high resolution measurements of the Ha line 
for the derivation of kinematic properties. The wave- 
length range also covers the [O l] A 6300, [N ll] A 6583 and 
[Sii] AA6717, 6731 A emission features. 

For the first SAMI commissioning run in July 2011, one 
field containing 13 galaxies was observed. The galaxies 
were selecte d from the 6dF Galaxy Survey (6dFGS; Jones 



et al. 2004) to have sufficient source density for SAMI, be" 
large enough to fill the spatial extent of the hexabundles 
(or, indeed, overfill them) and to have a surface bright- 
ness sufficient to give S/N>3 in the outermost ffbres in 
a 2-hour exposur e . The reader is referred to Section 4.2 
of [Groom et aT] ( |2012| ) for more details on the target 
selection. 

The observations are summarised in Table [H The data 
were observed using one of two strategies, either pointed 
or dithered, with the latter strategy being used to in- 



crease the spatial sampling and coverage of the observa- 
tions. The pointed data were observed on the 2nd July 
with a cloudless sky. The dithered data were observed 
on the 3rd July in less optimal conditions. For this pa- 
per the pointed and dithered data are kept separate for 
comparison. 

For the remainder of this paper, we concentrate on a 
single galaxy from the sample of 13 which were observed. 
This galaxy has 6dFGS identiffcation gl957222-550814 
and ESO identification ESO 185-G031 and is at a redshift 
of z=0.016. This corresponds to an approximate distance 
of 68 Mpc and an angular scale of \" equivalent to 340 
pc. An R-band image of the galaxy is shown in Figure 

m 

2.1. Imaging 

In addition to the SAMI data, ESO 185-G031 was 
observed with the Maryland-Magellan Tunable Filter 
(MMTF) and the f/2 camera of the Inamori-Magellan 
Areal Gamera and Spectrograph (IMAGS) camera on the 
6.5-meter Walter Baade Magellan telescope at Las Gam- 
panas Observatory on 22nd August 2011. The galaxy 
was observed with a Bessel R-band filter for three one- 
minute exposures. The MMTF order-blocking filter was 
then used at a central wavelength of 6660 A and a band- 
pass of 260 A. Four ten-minute exposures were taken 
in this configuration, yielding a total of 40 minutes on 
source. 

The combined three-minute R-band image of the 
galaxy is shown in Figure [l] The large red circle shows 
the position of the SAMI hexabundle for the pointed 
data. The dithers are not shown as they were small (of 
the order of 1/2 core diameter, or 0^8) and there were 
nine individual pointings observed. The smaller green 
circle, with a radius of \" shows the position in the galaxy 
of the supernova SN1998bo, a type Ic supernova seen in 
1998. T he d ata reduction and analysis are discussed in 
Sections 3.2 and 4.3 respectively. 



3. DATA REDUCTION 



data reduction package ( 


Groom, Saunders & Heald 2004 


Sharp & Birchall 20101), a robust package already used 



team to process SAMI data. It performs bias subtrac- 
tion and ffat-ffelding before extracting the individual fi- 
bre spectra. The extracted spectra are then wavelength 
calibrated, sky subtracted and corrected for fibre-to-fibre 
throughput variations, before being reformatted into a 
row-stacked spectrum (RSS) frame consisting of all 819 
extracted SAMI spectra (61 spectra for each of the 13 
hexabundles plus 26 separate sky spectra). 

For this data set, a spectrophotometric standard star, 
LTT6248, was observed on 2nd July. The star was ob- 
served in one hexabundle only and was reduced in exactly 
the same manner as the data frames. A "total" measured 
star spectrum is found by summing the extracted SAMI 
spectra. The measured star spectrum is compared to 
the tabulated fiux values for the star and a sensitivity 
function found for the SAMI spectra. This is normalised 
and applied to the extracted spectra, yielding relative 
fiux-calibrated RSS frames for each SAMI observation. 
The final step in the reduction is to scale each frame by 
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Figure 2. Left-hand panel: Integrated images of the sequence of dithered data cubes observed on 3rd July 2011. The images were formed 
by collapsing each data cube along the wavelength axis from 6300 A-7300 A. The shift from one observation to the next is clearly seen. To 
combine these data cubes it is necessary to calculate the positional offsets between each observation. Right-hand panel: The contours show 
the 2d Gaussian fit to the cross-correlation of each of the images in the left-hand panel with the fiducial image (in this case this is the first 
image in the sequence of dithers, so the first contour plot shown here is an auto-correlation). The postional offsets between the images are 
calculated from these fits. 



the exposure time to be in units of counts/second, thus 
making all frames directly comparable. 

The three pointed frames observed on 2nd July were 
then combined yielding a single RSS frame with a 2 hour 
exposure time and containing all thirteen galaxies in the 
field. 

3.1. The Dithered Data 

The SAMI hexabundles consist of a circular pattern of 
circular optical fibres. By its nature, the fibre coverage 
is not contiguous and, although the SAMI hexabundles 
have a high filling factor (greater than 75%), there are 
inevitably some gaps between the fibres. In addition, 
the circular fibres are not stacked in close-packed config- 
uration, but in a circular pattern, meaning the sampling 
also forms a non-uniform grid. These details present cer- 
tain challenges when working with dithered data. The 
non-contiguous, non-uniform coverage makes it difficult 
to simply interpolate over the holes in the sampling. We 
experimented with different ways to process the dithered 
data and our most consistent technique is described be- 
low. 

For the dithered data RSS frames were produced as 
usual using 2dfdr but before combining we first recon- 
struct entire data cubes for ESO 185-G031 from each 
SAMI observation. The SAMI data cubes used in this 
paper are created using a supersampled position grid at 
each wavelength slice. The sampling of the new grid is 
at the level of 1/lOth of a core diameter (0^16). The 
position of the centre of each core is defined to be the 
pixel closest to the tabulated position of that core, and 
so is determined on the new grid to an accuracy of 1/lOth 
of a core. A core is then defined to lie at that position 
in all wavelength slices. The core is "drawn" as an ap- 
proximate circle and filled with the spectrum recorded 



for that core. In this way each core is represented by a 
pixellated circle on the spatial grid. This represents the 
spatial extent of each sample without interpolating the 
data at this early stage. 

Once all data cubes are in hand they must be com- 
bined, accurately accounting for the positional offsets be- 
tween the cubes. To calculate the offsets between each 
cube in pixel space a cross-correlation method was used. 
First, each cube was collapsed over a large wavelength 
range (6300 -7300 A) to create an integrated image for 
each pointing. These images were then cross-correlated 
against a fiducial image, in this case simply the first 
SAMI image in our observing sequence. The resulting 
correlation function was then fit with a two-dimensional 
Gaussian to find the position of peak and therefore the 
offset from the peak position in the fiducial image. This 
allowed the data cubes to be accurately combined at the 
pixel level. The SAMI images and fitted cross-correlation 
functions are shown in the left-hand and right-hand pan- 
els of Figure [2] respectively. 

The cubes were combined by shifting and median 
stacking them using the offsets in pixel space from the 
cross-correlation method. This yielded a single combined 
dithered data cube. Note that the depth and therefore 
S/N of the dithered observation will vary across the field 
of view, again because of the non- contiguous sampling 
inherent in the hexabundles. 

For the analysis presented in this paper the combined 
cubes were then smoothed spatially with a Gaussian ker- 
nel with a 5-pixel sigma and spatially re-binned by a fac- 
tor of both two and four, yielding four final data cubes, 
two finely sampled (red and blue) and two more coarsely 
sampled (red and blue). The more coarsely sampled 
cubes are necessary for sufficient S/N per spatial pixel 
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in the blue and unless otherwise stated these are used 
for all further analysis. 

The effects of differential atmospheric refraction have 
not been corrected for in the f inal data cubes. However, 
following the formulation in Filippenko ( 1982 ) the effect 
is found to be small and does not impact the results of 
this paper. Within the individual line ratios the wave- 
length ranges are very short and the effect is negligible. 
Over the entire wavelength range considered when com- 
paring line ratios with one another (i.e. 4900 -6900 A) 
the relative shift is 0.35". This is half a pixel in the 
coarsely binned SAMI cubes (upon which this analysis 
was performed) and will therefore not affect our results. 

3.2. Imaging 

The MMTF data reduction pipeline was used to reduce 
the imaging data from IMMTF/MACS. Bias subtraction 
and flat-fielding were performed using Image Reduction 
and Analysis Facility (iraf) scripts and cosmic rays were 
identified and reject ed using the LACOSMIC routine of 
van Dokkum ( 2001| ). A median sky value for the entire 
image was calculated and subtracted. The astrometry for 
each image was found using the standard distortion map 
for IMACS, with some small adjustments. The images 
were registered using the USNO catalogue. The same 
procedure was applied to both the line and continuum 
images. 

4. RESULTS 

Once the SAMI data were reduced the analysis was 
performed using custom-written PYTHON scripts. The 
method was the same for both the pointed and dithered 
data. For each spectrum seven individual emission lines 
were each fit with a Gaussian function. These were 
the H/3, [Oiii] A 5007, [Oi] A 6300, Ha, [Nil] A 6583, 
[Sll] A 6717 and [Sll] A 6731 A lines. Where it was de- 
tected the underlying stellar H/3 absorption was also fit 
with a Gaussian and removed before fitting the emission 
line. Line strengths and kinematic properties were de- 
rived from the fits to each line. Six line strength maps 
derived from the dithered data are shown in Figure [3j 
Note that to produce the [S ll] map the line strengths of 
the two doublet lines were summed. From these six line 
maps it is possible to also derive four common emission 
line ratios (in logarithmic units), these being [Olll]/H/3, 
[OlJ/Ha, [Nil] /Ha and [Sll]/Ha. The spatial maps of 
these four ratios are shown in Figures |5] and [6j for the 
pointed and dithered data respectively. The"Nll]/Ha 
line ratio map is used to define spatial regions in the 
galaxy which are dominated by the disk and the galactic 
wind, using a threshold of [Nil] /Ha > —0.35 to isolate 
regions of the galaxy showing only ionisation due to star 
formation. Resulting integrated spectra are shown in 
Figure |4] 

4.1. Ionisation Diagnostic Diagrams 



First propos ed by Baldwin, Phillips fc Terlevich ( 1981 ) 
and revised byVeilleux & Osterbrock (l987'), Ionisation 
Diagnostic Diagrams (IDDs) are used to classify the dom- 
inant mechanism of gas ionisation in galaxies using the 
ratios of commonly observed optical emission lines. The 
most frequently used of these diagrams are those three 



Kewley et al. (2001) used stellar population synthesis 
models as well as photoionisation models to define a the- 
oretical maximum curve for star- forming objects on these 
three IDDs. Spectra of objects above this line (shown as 
the solid blue curves in our IDDs in Figures |7|and[8| 
are not dominated by star formation, but by a aifferent 
source of ionisation, for example AGN or LINER emis- 
sion. 



For the [Nil] IDD, [Kauffmann et al. | ( 2003D derived 
an alternative dividing line, based on empirical evidence 



from a large sample of galaxies. The [Kauffmann et 



al. 1(2003) curve (shown as^he dashed 
[IN llj IDDs in Figures 



our 



star- forming objects (below tne Kau"5mann et al. (2003) 



and \S\j 



lack curve m 
separates purely 



using the [O ill 
[Nil] /Ha and ' 



/H/3 line ratio versus each of the [O l]/Ha, 
S ll]/Ha line ratios. 



curve) from composite objects containin g both star for- 
mation and AGN/LINER-like activity (| Kewley et al.J 
"2006). Since we wish to separate the signature of star 
formation from any hint of other ionisation mechanisms 
(even if this ot her mechanism is not dominant) we adopt 
the ,Kauffmann et al. ( 2003[) line as the boundary be- 
tween star formation and AGN/LINER-like emission in 
the [Nil] IDD. 

4.1.1. The Pointed Data 

The IDDs for the pointed observations of ESO 185- 
G031 are shown in Figure [T] The green points corre- 
spond to individual spaxels for which both line ratios in 
question could be derived. (For some spaxels the S/N 
in the [O ill] and H/3 emission lines was too low to get a 
good fit. Since this ratio is required for all three IDDs, 
this is clearly the limiting factor here. See the final panel 
in Figure [5]) 

For the pointed data summed aperture spectra are also 
considered and plotted as yellow stars in Figure (T] The 
apertures are defined starting from the inner fibre and 
moving outwards, accounting for whole spaxels only, and 
roughly following the distribution of the [O lll]/H/3 emis- 
sion line ratio (see the final panel in Figure [sj where 
the numbers indicate spaxels summed to create the six 
aperture and annular spectra). The aperture spectra are 
found by summing the spectra within the aperture and 
the line ratios are derived as usual. As can be seen in 
Figure [t] all integrated aperture spectra (yellow stars) 
lie very close to or on top of one another. This shows 
that the light from the star forming disk is so dominant 
that integrated spectra over a fixed aperture, even if that 
aperture was the size of the entire SAMI hexabundle, 
would not pick up the non-star-forming signature we see 
off the plane of the disk. 

The red diamonds in Figure [7| show the same apertures 
but now considered as annuli, i.e. this time each annulus 
is plotted individually instead of summing over all light 
within the aperture in question. This time a clear trend 
is seen, with the central annuli moving along a track 
through the star-forming region of the diagram towards 
the limiting curves and the outer two spectra clearly ex- 
hibiting properties indicative of a non-star-formation ion- 
isation source. 

Again, it is clear that with a single fibre or a single large 
aperture these characteristics would have been missed. 
It is the spatial sampling of our data that allows us to 
discover the true nature of this galaxy. 

4.1.2. The Dithered Data 
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Figure 3. Spatial maps of the emission line strengths derived from the dithered data where the [S ll] map is a sum of the two lines i 
doublet. The disk is clearly seen in emission in all lines. The colour bar is in units of logio(counts/s). The images are 24^^on a side. 
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Figure 4. Left: two integrated SAMI spectra, from the disk and wind regions of ESO 185-G031, covering the blue wavelength range 
from 4900 -5100 A. The and [O ill] emission lines are clearly seen in both spectra, but with lower S/N in the wind region. Right: two 
integrated SAMI spectra, from the disk and wind regions of ESO 185-G031, covering the red wavelength range from 6400 -6850 A. The 
[Ol], [Nil], Ha and [S ll] emission lines are clearly seen in both spectra, but with lower S/N in the wind region, where the sky residuals are 
also more prominent. 



The IDDs for the dithered d ata are shown in the top 
part of Figure [sj A s weh as the |Kewley et al. (2001 ) and 
Kauffmann et al.] (12003) Unes discussed above , we also 
present two tidu^aHme s following the work of |Sharp~fcl 
Bland-Hawthorn] (|2010[ ). These are shown as solid red 
lines in all three IDDs. The shallower line is a fit to 
the track followed by points in NGC1482 with increasing 
height above the plane of the disk of the galaxy. Since 
NGC1482 is a canonical starburst-driven galactic wind 
( [Veilleux fc Rupke ||2002[ ), this line is taken to represent 



the trend expected in objects of that kind. The points lie 
in this part of the diagram because the ionisation source 
is shock excitation caused by the starburst-driven wind 
in the galaxy. Likewise, the steeper line is a fit to the 
track followe d by points in NG C1365, a canonical AGN- 
driven wind ( Veron et al. |198Q ) and is taken to represent 
the trend expected in objects of that kind. The solid 
black line is the bisector of the two lines and is used to 
partition points into one or other category, where they 
are colour-coded: blue for AGN and red for starburst- 
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Figure 5. Spatial maps of the emission line ratios derived from 
the pointed data. Blank spaxels correspond to spectra where the 
S/N in the emission lines in question was not high enough to obtain 
a good fit. These maps are used to construct the IDDs in Figure 
Since the [O lll]/H/3 ratio (final panel) is needed for all three ID! 
incompleteness in this map clearly has the greatest impact. The 
final panel shows the spaxels numb ered in grou^ used to extract 

maps a 
! plane of the disk of 
the galaxy, indicating the presence of a galactic wind. The images 
are 24^^on a side. 



final panel shows the spaxels numb ered in groups used to 
aperture spectra (see Section [4.1.1| and Figure m). In all 
trend towards higher line ratios is seen off the plane of the 



[Nlll/H, 



[Sll]/Ha 




Figure 6. Spatial maps of the emission line ratios derived from 
the dithered data. Again, the blank spaxels are those without 
spatial coverage or sufficient S/N to fit the necessary lines. These 
are used to construct the IDDs in Figure [s] The trend towards 
higher line ratios off the plane of the disk is clearly seen, indicating 
the presence of a galactic wind. The images are 24^^on a side. 

driven shocks. 

The bottom part of Figure [S] shows the colour-coded 
points re-mapped onto their spatial position in the 
galaxy. This helps to determine the spatial origin of the 
different ionisation mechanisms which are potentially in 
play. In all three IDD maps the disk of the galaxy is 
dominated, as expected, by star formation. However, the 
picture becomes more complex as we move off the plane 
of the disk. In all three IDDs an indication is given that a 
different, harder, ionisation source is at work away from 
the disk. The [N ll] IDD and map are ambiguous in that 
they show evidence for both AGN excitation and shock 
excitation in the south-west and north-east parts of the 
galaxy. 



The [S ll] diagram is less ambiguous pointing towards 
both wind regions being dominated by ionisation from 
starburst driven shocks. The [Ol] diagram also hints 
more strongly towards shock ionisation over AGN ioni- 
sation. For these reasons the starburst-driven wind sce- 
nario is favoured, but it is not possible to entirely rule 
out a wind driven by an obscured AGN on the basis of 
the IDDs only. 

The star formation rate (SFR) for ESQ 185-G031 is 
derived from IRAS 6Q/im (Fen) and 100/im (Fioo) flux 

([2002). The mea- 
1.66 Jy, yield- 
1.92 X 



measurements after |Kewley et al. 
sured values are F6o=0.52Jy and r loo 
ing a value for far-infrared luminosity of L(FIR) 



10 ergs . This luminosity represents that over the 
range 40-120/im. The ratio of L(FIR) to L(IR), the 
latter being the tot al IR dust em ission, is given by 
L(IR)=1.75L(FIR) ( [Calzetti et al.||2000| . This then 
gives a SFR of I.TMqYT-' for ESQ 185-G031. 

The use of this conversion between L(FIR) and L(IR), 
and of L(IR) to calculate star formation rate, assumes 
that all dust emissi on is due to heati ng from the young 
stellar population (Calzetti et al. |20QQ[ ). If this is not the 
case, and some heating is from other sources such as an 
older population of stars, then this will overestimate the 
SFR. However, assuming the converse, that L(FIR) is a 
more accurate representation of the bolometric luminos- 
ity, will certainly underestimate the SFR. With this in 
mind we adopt 1.7M0yr~^ as the correct SFR for this 
galaxy. 

A second, independent, measurement of SFR can 
be made using ultraviolet emission. GALEX observa- 
tions of this galaxy give mFuv= 19-05 mag wi th E(B- 
V) =0. 054 mag. Using the procedure in Lopez-Sanchez 
( |2010[ ) yields a value for SFR of 0.07 M© yr-\ A poten- 
tial source of this discrepancy is the estimate of obscu- 
ration from the UV data. Given the brightness of this 
galaxy in FIR the true obscuration is likely to be larger, 
meaning the SFR estimate from the ultraviolet should 
be taken as a lower limit. The difference between these 
two indicators could also indicate the importance played 
by the star formation hist ory of the galaxy in determin- 
ing an accurate SFR (e.g. Koribalski fc Lopez-Sanchez 

However, individual 



2009 Lopez-Sanchez et al. 



estimates of SFR are still likely to be discrepant simply 
because the calibration of different SFR indicators can 
give estimat es which vary by up to a factor of two (see 
for example 
From the 



Hopkins et al.|( 2003D). 
MAGS/MMTJ 



imaging it is estimated that 
the star formation is occurring in the two prominent 
Ha "eyes" covering an area roughly 3kpc in diame- 
ter. This implies a star formation surface density of 
about O.24M0yr-ikpc-^ This is more than sufficient 



to launch a bipolar wind ( Heckmaii|[2002 ) . In addition 
this g alaxy does not appear in the SUM SS (Mauch et al. 
|2003[ ) or AT20G ( iMurphy et al. |2010[) catalogues imply- 
mg that it is unlikely to host an AGN and the observed 
large-scale wind is therefore most likely driven by star 
formation. 

An estimate of the metallicity of the star-forming disk 
was derived using only those spaxels showing a signature 
of star formation. We derived an integrated spectrum of 
the disk by masking the data cube based on the [N ll]/Ha 
line ratio, rejecting spaxels with [Nil] /Ha > —0.35 and 
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summing the resulting cube. We calculated values of 
[Oiii]/H/3=0.577±0.048 and [N iil/Ha=0.349±0.018 and 
used the relations in (Pettini fc Pagel^ 2QQ4) to empiri- 
cally determine the gas-phase oxygen aWndance of this 
galaxy. We find a consistent value of 12+log(0/H)=8.64 
and 8.66 dex using the N2 and the O3N2 parameters, 
respectively. 

It is important to note that these oxygen abundances 
correspo nd to values calibrated by follo wing the direct 
method (jOsterbrock & Ferland 1 12006|). To compare 



Figure 7. IDDs derived from the pointed observations of ESO 185-G031. The 'Kewley et al~| ( |200 1^ hmiting starburts curves are shown 
for ah three IDDs (sohd blue curve) with the Kauffmann et al. (2003) delimiting line also shown tor the [Nil] IDD (dashed black curve). 
The yellow stars show the position of the galaxy for incrementally summed apertures, all of which indicate a system entirely dominated 
by star formation. The red diamonds show the same apertures treated as incremental annuli, with the outer two annuli labeled (5 and 6). 
This time a trend is seen such that the outer parts of the galaxy (annuli 5 and 6) show ionisation dominated by shocks from the newly 
discovered galactic wind. This clearly demonstrates the need for spatially-reolved spectroscopy to identify such systems in the future. A 
representative error bar is shown in green in the top-right corner of each plot. 

the Ha emission line is high enough to fit in all spaxels. 
The dithered observations reveal the trends in the emis- 
sion line kinematic properties much more convincingly. 
The top left panel of Figure[TO shows the Ra velocity field 
extracted from the dithered data. The structure in the 
map is very clear with the kink in the south of the disk 
readily apparent. The bottom right panel shows the de- 
rived FWHM of the Ha emission line. The trend seen in 
the pointed data is reproduced, with significantly broader 
emission lines seen off the plane of the disk in both di- 
rections. The regions of the data with disturbed velocity 
and broadened emission lines correspond to those regions 
displaying non-star-formation driven ionisation proper- 
ties. 

The shape of the emission line changes off the plane 
of the disk as well, hinting at more than one velocity 
component in the gas (e.g. rotating disk and wind out- 
fiow) . To illustrate this we extract spectra from four 2x2 
spaxel apertures, as shown in the bottom right panel of 
Figure [To] The corresponding spectra (zoomed in on the 
Ha emission line) and Gaussian fits are shown in the 
four panels of Figure [TT] Apertures 2 and 3 are in the 
plane of the disk and show significantly narrower line pro- 
files (FWHM=77±7 and 59±3kms-^ respectively) than 
those seen in apertures 1 and 4 (FWHM=119±13 and 
130±16kms~^ respectively), which lie in the "wind" re- 
gions. This line-broadening supports the idea that we 
are observing two separate kinematic components super- 
imposed, the rotating, star-forming disk and the wind 
outfiow. 

To further examine the nature of the observed velocity 
field a rotating disk model was used to fit the obser- 
vations. Before the fit was performed the Ha velocity 
map was masked using the [Nil] /Ha map to distinguish 
the parts of the map corresponding to star-forming disk. 
A threshold of [Nil] /Ha > —0.35 effectively masks the 



with oxygen abundances calculated using photoioniza- 
tion models, these values should be increased by 0.3- 
0.4 dex, i.e. the gas-p hase metallicity of th is galaxy is 
12+log(O/H)-9.0 (see [Lopez-Sanchez et al ||2Q12[ for a 
review of this problem). 

4.2. Kinematics 

4.2.1. The Pointed Data 

Gas kinematics for the pointed data were derived using 
a Gaussian fit to the Ha emission line in each spaxel. 
Figure [9] shows the results, with a quasi-continuum map 
(summed from 6300 - 7300 A) in the first panel, the Ha 
line strength map in the second and the velocity and 
FWHM maps in the third and fourth respectively. 

Close examination of the velocity field reveals some 
asymmetry, most notably a kink in the redshifted side 
of the velocity map, south of the plane of the disk. The 
sigma map also shows a trend as one moves off the disk, 
with the emission line broadening in the gas off the plane 
of the disk. Broadening in this manner can be due to two 
distinct velocity components, such as a rotating disk and 
an outfiow, superimposed along the line of sight, thus 
supporting the existence of a wind in this galaxy. 

4.2.2. The Dithered Data 

For the dithered observations the Ha kinematics are 
derived using the finely binned data cubes as the S/N in 
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Figure 8. Top panels: IDDs derived from the dithered observati ons of ESQ 185-G031. Ag ain, the [Kewley et a l."|(|2001) limiting starburts 
curves are shown for all three IDDs (solid blue cu rve) with the 'Kauffmann et al. | ([2003 ') delimitmg Ime also shown for the [N ll] IDD 
(dashed black curve). The colour coding is that of [Sharp Bl and-Hawthorn (| 2UlO| ), where green denotes spatial points dominated by 
star formation, blue denotes a dominant AGN-like emission and red denotes shock emission characteristic of a starburst driven wind. A 
representative error bar is shown in black in the top-left corner of each plot. Bottom panels: The colour-coded points are re-mapped to their 
spatial positions in the galaxy. All three maps show a star-formation-dominated disk with a stron ger ion isation mechanism dominating off 
the disk in both directions. This is strong evidence for a galactic wind in this galaxy. See Section |4.1.2[ for more details. 



parts of the galaxy showing non-star- forming ionisation. 

We use a simple disk model to fit the projected Ha 
velocities. Our model uses a n arctangent rotation curve 
( St aveley- Smith et al. |p!99Q ), given by 



V{r) 



asym , 

— - — arctan 

TT 



(1) 



where Ksym is the asymptotic circular velocity, and Vd 
is the kinematic scale radius. The position angle and 
inclination determine the 3D spatial orientation of the 
disk. The model also includes the position of the cen- 
tre of the disk within the velocity map, and a system- 
atic velocity offset (which corrects for any bias in the 
systemic redshift) as free parameters. The model pro- 
duces a velocity map analogous to the observed map. 
A Levenberg-Marquardt minimisation routine optimises 
the model's seven free parameters to give the best fitting 
model velocity map. 

The results are shown in Figure [lOj with the disk model 
displayed in the top right panel and the difference be- 



tween the observed velocity map and the model in the 
bottom left panel. The residual map shows the region in 
the south-west of the observations contains gas not rotat- 
ing with the disk. We interpret this as material flowing 
away from us in an outflow. The outflow may in fact 
be bipolar as is hinted in the Ha FWHM map but we 
do not have enough spatial coverage to the north of the 
galaxy to confirm this. Of note as well is the sharpness of 
the velocity transition between the rotating disk and the 
extra-planar material. This is thought to be due to the 
inclination of the disk, which is derived from the velocity 
model to be 63 degrees. With the disk inclined thus to 
the line of sight the dusty edge of the disk is seen to the 
south-west of the disk (see Figure [T]) and probably corre- 
sponds to the region of low FWHMin the Ha kinematics 
(see panel 4 of Figure 10). 

The fit to the Ha kinematics is taken to represent the 
genuine rotation of the gas in the galaxy. Velocity maps 
are derived from all other emission lines and the disk 
model is subtracted from each in turn to examine de- 
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Figure 9. The Ho; kinematics for ESO 185-G031 derived from 
the pointed observations. A kink is seen in the south-west of the 
velocity map and a trend towards broader hne widths off the plane 
of the disk is seen in the FWHM map. These are both indicative 
of a complex kinematic structure comprised of more than one kine- 
matically distinct component, supporting the presence of a wind 
in this galaxy. 
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Figure 10. The Ra velocity map for ESO 185-G031 (top left) 
shown with the best fit disk model (top right; see Section |4.2.2| 
for details). The residual map, found by subtracting the rnodel 
from the observed velocity map (bottom left), shows a very strong 
deviation from the rotating disk model in the region of the galactic 
wind. This region is co-spatial with the broadened emission lines 
seen in the FWHM map (bottom right). The squares overlaid 
on the FWHM map correspond to the four apertures for which 
individual Ha spectra are displayed in Figure pT] 

viations from the rotating disk. For the red emission 
hnes, [Oi]A6300, [Nn]A6583 and [Sh]A6717A the finely 
binned cubes are used. The results are shown in Figure 
12 In all three velocity maps the signature of the wind 
is clear. For the blue emission lines the more coarsely 
binned cubes are use to derive the kinematic properties 
(with the disk model binned accordingly before subtract- 
ing). The results are shown in FigureflS] Again, the wind 
signature is seen in the south-west oitne residual maps. 

4.3. Imaging 

For the IMACS/MMTF image analysis, each of the 
individual frames were aligned on stars, scaled to a con- 
stant sky value, and the stack mean was used to form 




6660 6665 6670 6675 6660 6665 6670 6675 6680 

Wavelength (A) 

Figure 11. Four spectra extracted from 2x2 spaxel apertures. 
The extracted spectra are blue, with the single Gaussian fits shown 
in green. Apertures 2 and 3 in the disk show narrower Gaussian 
profiles, than apertures 1 and 4 (in the "wind" part of the galaxy) 
which show significantly broader profiles. This indicates that the 
emission in this region may come from two kinematically distinct 
components - the disk of the galaxy and the wind. 
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Figure 12. The left-hand column shows the velocity maps derived 
for each of the [Nil] A 6583, [Sii]A6717 and [Oi]A6300A emission 
lines. The right-hand column shows the residuals found when the 
disk fit to the Ho; velocity map is subtracted from each of the maps. 
The wind signature is clearly seen in each of the residual maps. 
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Figure 13. The left-hand column shows the velocity maps derived 
for the and [O ill] A 5007 A emission lines. The right-hand col- 
umn shows the residuals found when the disk fit to the Ho; velocity 
map is subtracted from each of the maps. The wind signature is 
clearly seen in each of the residual maps. 

the science image except where the brightest pixel in the 
stack deviated strongly from the rest so was replaced by 
the median before taking the stack mean. Seven cycles 
of Lucy-Richardson deconvolution were applied to the 
line image to match its stellar PSF to that of the contin- 
uum image. The continuum was scaled to produce zero 
subtraction residuals within the galaxy; a scale change 
of 2.7% between line and continuum was also applied. 
Foreground stars are still evident because of the MMTF 
passband used here. In Fig. [T4j two Ha-bright "eyes" 
straddle the galaxy center in the disk, while many faint 
filaments extend above the disk. This extended emission 
strongly supports the discovery of a wind in this galaxy. 

5. CONCLUSIONS 

We have presented SAMI IFU observations of a single 
galaxy observed in July 2011. This galaxy, ESO 185- 
G031, was part of the commissioning sample of 13 and 
showed interesting structure in its kinematic and optical 
line ratio maps. All four of the measurable line ratios 
show a significant trend to higher values off the plane of 
the disk. Using IDDs to further investigate this reveals 
that emission in the disk is dominated by ionisation from 
star formation, but the regions off the plane of the disk 
exhibit ionisation by some other mechanism. We pro- 
pose that this other mechanism is shock excitation from 
a starburst-driven wind, based on comparison with the 
IDDs for NGC1482, a canonical star burst driven wind 
( [Sharp fc Bland-Hawthorn | |2Q1Q|). 

The galaxy also shows deviations from a simple rotat- 
ing disk in its spatially resolved velocity field, implying 
that more than one kinematically distinct processes are 
at play. This is borne out when comparing the measured 
velocity field to a disk model fit to the part of the velocity 
field corresponding to the star forming disk (as defined 
by the [Nil] /Ha emission line ratio). A clear velocity 




Figure 14. A 31^^ x M" region extracted from the MMTF im- 
age, where north points upwards and east to left. This continuum 
subtracted emission-line-only image is a composite, with inverted 
gray scale to show the outer parts of the galaxy smoothed to 0(^8 
FWHM and positive gray scale showing the disk at 0(^6. Numerous 
faint Hq;+[Nii] filaments extend up to 11" (3.7 kpc) from the disk. 



residual is seen which is co-spatial with the region of el- 
evated line ratios, supporting to the wind hypothesis for 
this galaxy. 

Previous studies have identified winds_through__gas 
kinematics using long-slit spectroscopy (Axon & Taylor 



1978) or through variations in gas excitation using tun 
able filters (Veilleux & Rupke 2002). Here we demon- 
strate that integral field spectroscopy is able to identify 
galactic winds through a combination of gas excitation 
and kinematics. Thus we anticipate that a large survey 
using hexabundles will enable the first objective survey 
of the prevalence of outflows in galaxies. 

More generally, SAMI will be very effective at demon- 
strating departures from symmetry over a wide range of 
galaxy types, both in the stars and in the ionised gas. We 
anticipate that all or most galaxies will revea l such be- 
haviou r, as evidenced by the recent SAURON Jde Zeeuw 
jet al. (2002) and ATLAS-3D (|Cappellar i et al. |2Ulip sur- 
veys. Future samples will be so large that we can study 
this behavi our as a function of environment across the 
hierarchy ( Lawrence et al. ||2012 ). Just how galaxies are 
shaped by their environments is one of the great un- 
knowns of modern astronomy. We fully anticipate that 
hexabundle survey spectrographs will have a major im- 
pact in answering this question in the years ahead. 
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